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Abstract 

We study hydrogen-like atoms in Af = 1 supersymmetric quantum electrody- 
namics with an electronic and a muonic family. These atoms are bound states 
of an anti-muon and an electron or their superpartners. The exchange of a 
photino converts different bound states into each other. We determine the en- 
ergy eigenstates and calculate the spectrum to fourth order in the fine structure 
constant. A difference between these perfect atoms and non-supersymmetric 
ones is the absence of hyperfine structure. We organize the eigenstates into 
super multiplets of the underlying symmetry algebra. 
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1 Introduction 



Supersymmetry is often invoked to resolve a number of theoretical difficulties with the stan- 
dard model of particle physics. The symmetry can control quantum corrections to the Higgs 
mass, thus providing a solution to the hierarchy problem. The symmetry suggests the strong 
force, the weak force, and electrogmagnetism are unified at high energy scales. Moreover, the 
lightest supersymmetric partner is a candidate for dark matter. Despite these theoretical ad- 
vantages, we see no direct evidence for supersymmetry at low energies; supersymmetry must 
be broken, and most studies of supersymmetry are devoted to investigating methods for and 
consequences of the breaking. In this paper, we take a different tack and look at the energy 
spectrum of an anti-muon electron bound state in a theory with unbroken supersymmetry, 
supersymmetric quantum electrodynamics (SQED). 

We were initially inspired to write this paper by work on gauge/gravity duality. A 
gauge/gravity duality is a map between a field theory and a string theory. The duality 
is useful because when the field theory is strongly interacting, the string theory is weakly 
interacting and vice versa. Both the field theory and the string theory are typically su- 
persymmetric. Often one is faced with the following awkward situation: A calculation on 
the gravity side has revealed some property of the strongly interacting field theory, and the 
corresponding property of the field theory at weak coupling has not yet been studied. One 
prime example of such a situation was the computation of the viscosity of M = 4 SU(iV) 
super Yang-Mills theory at strong coupling in ref. pp. Only five years later was the viscosity 
calculated in the perturbative limit [2]. In the case of supersymmetric atoms, ref. [3JH] stud- 
ied hydrogenic bound states at strong coupling in A" = 4 SU(iV) super Yang-Mills modified 
by the addition of two massive A~ = 2 hypermultipletsQ No corresponding study at the time 
had been made of such bound states at weak coupling. Moreover, the interesting observation 
was made that these bound states exhibited no hyperfine structure jl]. The absence of such 
structure is an almost trivial consequence of Af = 2 supersymmetry [5], but in this paper 
we will see that hydrogenic atoms of Af = 1 SQED also lack hyperfine structure. It should, 
however, be emphasized that the energy spectra of A/" = 1 and Af = 2 SQED remain notice- 
ably different. The energy levels of Af = 2 hydrogenic atoms are independent of the spin 
of both the electron and the proton, while fine structure effects remain evident for Af = 1 
atoms. 

A second motivation for this paper is pure intellectual curiosity. Although the Is state of 
supersymmetric positronium was considered almost thirty years ago [B] , no one to our knowl- 
edge has studied anti-muon electron bound states in SQED. The way in which the bound 
states organize themselves into supermultiplets is surprising and intricate. Similar to what 
happens for supersymmetric positronium, both degenerate and second order perturbation 
theory contribute at the same order in the fine structure constant a. 

We hope that these super atoms may be useful in particle physics, perhaps as a candidate 
for dark matter, perhaps in a hidden sector, perhaps for neutrino physics. 

The paper is organized as follows. In section |2l we review SQED. In section [31 we present 
the relevant scattering amplitudes necessary for computing the energy spectrum to order 
a 4 . In section HI we reduce the energy spectrum computation from field theory to time 

1 Given the strongly interacting nature of the bound states, heavy- light or hybrid meson is perhaps more 
appropriate terminology. 
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independent perturbation theory in quantum mechanics. Section [5] contains detailed results 
for the hydrogenic states and their energies. Section [H] contains some discussion. First, 
however, we summarize our results. 

1.1 Results 

The energy spectrum of the hydrogen atom is usually described order by order in the fine 
structure constant a. The rest mass of the atom is M + m where M is the mass of the proton 
and m the mass of the electron. The binding energy is of order a 2 fi where \x = MmjiM +m) 
is the reduced mass and we work in units where the speed of light c = 1. Fine structure 
effects are of order a 4 /z and involve relativistic corrections along with spin-orbit couplings 
of the electron's spin to its orbital angular momentum. Hyperfine structure is of order 
a A nm/M and involves spin-spin coupling of the electron and proton. There are higher order 
corrections, for example the Lamb shift at order a 5 , but in this paper we work only to order 
a 4 . 

The proton is a composite object in the real world, and its compositeness has subtle 
effects on the hydrogen spectrum that do not interest us for the purposes of this paper. 
Thus, we replace the proton with a fundamental particle of positive charge, an anti-muon. 
Although in the case of the hydrogen atom M is much larger than m, the results we present 
are valid for arbitrary values of m and M. 

We consider SQED in 3+1 dimensions with four super charges. The electron and muon 
have super partners, the selectron and smuon. Because the electron and muon have both 
charge and a Dirac mass, they need to be Dirac fermions and as such will each have two 
complex scalar field super partners. In other words, there are two selectrons and two smuons. 
The super partner of the photon is a Majorana spinor, the photino. 

The existence of these super partners leads to some interesting effects. An electron anti- 
muon bound state can mutate into a selectron anti-smuon bound state and back through 
photino exchange. There are also fermionic bound states: an electron anti-smuon or a se- 
lectron anti-muon which can mutate into each other. Because of photino exchange, the 
eigenstates of the effective Hamiltonian describing our super atom are actually linear super- 
positions of these different types of bound states. 

The total angular momentum is a good quantum number for the bound states and is a 
useful organizing principle for the energy spectrum. Consider a hydrogenic wave function 
with principal quantum number n and orbital angular momentum I > 0. Let VJ be a 21 + 1 
dimensional representation of the SO (3) rotation group. A bound state of an electron and 
an anti-muon will transform as 



under 50(3). As there are two selectrons and two smuons, there are four fermionic bound 
states consisting either of an electron anti-smuon or selectron anti-muon. These bound states 
transform as 



Finally, there are four bosonic bound states consisting of a selectron and anti-smuon, all 
transforming as VJ. 



Vi <8> V l/2 g> V l/2 = Vi-x ®2Vi@ V l+1 
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Given supersymmetry, the energy spectrum must organize itself into super multiplets 
arising from the four super charges in 3+1 dimensions. If j is a total angular momentum 
quantum number, then for j > 0, a massive super multiplet consists of the four representa- 
tions of the rotation group IZj = Vj-1/2 © 2Vj © ^j+i/2- From the analysis in the previous 
paragraph, we see that to each hydrogenic wave function of principle quantum number n 
and orbital angular momentum I > 0, we can associate the four super multiplets 
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Our results for the first few atomic energy levels of supersymmetric hydrogen are given 
in Figure [TJ More generally, we find that the two multiplets TZi are degenerate in energy. 
The multiplet 1Zi + \/2 is higher in energy by an amount 
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2(1 + 1)(2/ + l)n 3 ' 
while the multiplet Hi-1/2 is lower in energy by an amount 
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In the case j = 0, the H-1/2 multiplet of course does not exist and the TZq multiplets do not 
contain the representation V-1/2' 
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The two TZq multiplets are degenerate in energy, and the 72-1/2 multiplet is higher in energy 
by an amount /ia 4 /2n 3 . 

There is also an overall shift in the energies at order a 4 , see ( 14. 6p . For the total energy 
of a state in a super multiplet IZj, we find 
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As is also true in QED, this expression can be written purely in terms of j. In other words, 
the super multiplet TZi + i/ 2 that comes from an \nl) state and the 1Z(i+\)-i/2 super multiplet 
that comes from an \n, I + 1) state are degenerate in energy. 

An important observation about these energy splittings (11. Ah and (11.51) is the absence 
of hyperfine structure; there is no energy splitting between two states of order a 4 /j,m/M in 
the case m <C M. This absence is an effect of supersymmetry. In pure QED, we find the 
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Figure 1: Spectrum. Splitting of the energy levels n = 1,2,3 with angular momentum Z = 0, 1, 2. 

2 

The labels on the lines indicate the super multiplets. The bare energies — ^fr are shifted due to a 
relativistic correction to the kinetic energy by 5E n [ as given in (|4.6p . and then split due to various 
interactions into two or three levels. 
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multiplets VJ © V i+1 = A l+1/2 (I = 0,1,2,...) and Vj_i © Vj = Ai- 1/2 (I = 1,2,...). In the 
limit m <C M, there is a fine structure splitting between Ai + \/2 and Ai-1/2, and a further 
hyperfine splitting between the V\ and Vi±\ in Ai±\/2- 

In cases with more supersymmetry, the corresponding multiplets are even larger and 
the splitting (11 -4p and (jl.5p will disappear as well [5]. For example in a theory with 8 
supercharges, a massive multiplet with j > 1 transforms as Vj_iffi4V^_i/ 2 ©6Vj©4Vj +1 / 2 ffiVj-+i 
under the rotation group. To get such a large multiplet, we need to combine the four 
multiplets in table (II. 3p . Similarly, in the case j = 0, we would need to combine the three 
multiplets in table (jl.fip . 

It is also instructive to compare our results for muonium to the splitting of the ground 
state energy of Af = 1 positronium found in [6j. The two computations differ in the respect 
that the latter one involves two particles of the same mass that moreover can annihilate. 
Setting M = m and n = 1 in (11. 7p . we find that the two levels at Z = differ in energy 
by AE = which is half the value for the splitting between the ortho and para states of 
positronium j6j. This difference is a consequence of the absence of the annihilation diagrams. 

2 SQED 

We write M = 1 SQED for two families of matter particles which we refer to as "electronic" 
and "muonic." The electron e~, its superpartners the selectrons e±, and their antiparticles 
e + and e±, are collectively represented by two chiral superfields <& e ± with U(l) charge ±e 
and mass m. Similarly, we write $ m ± for the muon (/z~,/! + ) and the smuons which 
are assigned mass M. The U(l) gauge superfield containing the photon 7 and the photino 
7 is denoted by V. Employing the superspace conventions of Wess and Bagger |7|, the 
Lagrangian has the form 

£sqed = \{ww\ e2 + ww\- e2 ) 

+ (<Pi + e 2eV <P e+ + <Pl_e- 2eV <P e _)\ e2§2 +m(^> e+ <P e _\ 02 + 

+ (< + e 2eV m+ + <_e- 2eV m „) \ e2§2 + M($ m+ $ m \ 2 + $l + $l_\- e2 ) ,(2.1) 

where the super fieldstrength is defined by W a = —^D 2 D a V and W& = —\D 2 DaV. After 
integrating out the auxiliary fields, the Lagrangian can be written as a kinetic term for the 
gauge fields 

£gau gc = —F^F^ + % -\Yd,\ , (2.2) 
a part that contains the electronic particles 

^electron = ^e(*7^ + ™) tfj e + 4 + (V 2 - m 2 )(j) e+ + <\>\_ (l? 2 - m 2 )<f) e _ 

+ V2ie(<j) e+ 4> e P^\ - (pl+XP+ipe - e _AP_V e + <pi_4> e P+\) , (2.3) 

an analogous one for the muons, £ mu0 n, which is obtained by replacing the labels e by m, 
and a part with contact interactions between the two families 

6 2 

Contact = -^{\4>e+\ 2 ~ \4>e-\ 2 + \<Pm+\ 2 ~ \4>m~\ 2 ) ■ (2.4) 
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For our notation and conventions, see App. |A] 



3 Scattering amplitudes 

In order to find the spectrum of bound states of a particle of the electronic family and an 
anti-particle of the muonic family in Sec. [5j we first compute the potential between any 
two of these particles. We deduce the potential energies from the non-relativistic limit of 
the tree-level scattering amplitudes which we compute from SQED Feynman rules. The 
amplitudes will allow us to calculate the bound state spectrum including all effects up to 
order a 4 in the fine structure constant a — f-. 

4-7T 

In the next subsection, we will explicate the derivation of the potential from the ampli- 
tudes for the scattering of an electron and an anti-muon. This scattering process is the only 
one that would exist for pure QED. The amplitudes and results for all other cases are listed 
in the two subsequent subsections. 



3.1 QED 

At tree-level the only diagram describing the scattering of an electron and an anti-muon 
involves the exchange of a photon: 

) = e 2 u e (q)^u e (p) A^ u (p - q) v m (p)>y u Vm(q) • (3.1) 

H + * • + /i + 

The in-going electron and anti-muon momenta are p and p' respectively. The outgoing 
momenta are q and q'. The spinors are 



u,( P ) = I , „.(„) = ^4f- I , (3.2) 



CO 



Imo 

We will work in the center of mass frame where 



p = (\fm 2 +p^,p) , q = (\/m 2 + q 2 , q) , (3.4) 

p' = (^/M 2 +p 2 , -p) , q' = (^M 2 + q 2 , -q) . (3.5) 
The first few terms in a non-relativistic expansion of the scattering amplitude are^ 



4Mme 2 . + \ 1 2 fl 1 \ 2 (p 2 - q 2 f 1 



2^ ^ \2m 2 Mm) 2^ ^ \2M 2 Mm, 
X 1 
+ 4Mm~^~ ^ 2<?e ' dm ~ AMm^~ £M e (p - <£) • <r m + . . . VmoU • 



2 We are grateful to Tomas Rube and Jay Wacker for pointing out a mistake in this and similar formulas 
in a previous version of the manuscript. Before we worked in Feynman gauge where it would have been 
necessary to compute a one- loop diagram to fix an ambiguity in the non-relativistic potential [8] . 
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The ellipses denote terms that are higher order in the space-like momenta p and q. The 
plane wave states in quantum field theory are normalized to the Lorentz invariant quantity: 

(p\q) = 2^m 2 + q 2 (2vr) 3 5 (3) (p - q) . (3.7) 

In non-relativistic quantum mechanics, in contrast, these plane wave states are typically 
normalized to (2ti) 3 S^(p — q). To take into account the different normalizations, we define 
the non-relativistic scattering amplitude 

M 



-Mnr = — . , 

4 [(m 2 + p 2 ) (m 2 + q 2 )(M 2 + p 2 ) (M 2 + q 2 ) 1 

Taking into account the change in normalization, we find that 
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(3.9) 



We have changed the spinor 77 of the anti-muon into a spinor £ as if it described a muon, 
77 = icr 2 £*. 

We would like to compare this scattering amplitude with the Born approximation result 
for a particle of position r and momentum p in non-relativistic quantum mechanics scattering 
off of a potential V(f,p). The Born approximation says that 
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(3.10) 



for plane wave initial and final states. We now Fourier transform the amplitude .Mnr with 
respect to p — q, keeping p as a variable. We find FT(.Mnr) = — V(r, p) and 
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(3.11) 



All terms are understood to be normal ordered, i.e. when p and L are replaced by operators 
then they do not act on the coordinate dependence of the potential. The result (13.111) is 
familiar up to subleading corrections in 1/M. The first term is the Coulomb attraction. The 
second term is the orbit-orbit, also referred to as the current-current, coupling. The third 
term is the Darwin term. The fourth term is the spin-orbit coupling of the electron. The fifth 
term is the spin orbit coupling of the muon. The last term is the hyperfine coupling between 
the spin of the electron and the spin of the muon. For a hydrogenic orbital, the expectation 
values of (1/r) and (p) scale as ct/x Thus, this non-relativistic expansion of the effective 
potential is also an expansion in the fine structure constant. The Coulomb interaction is of 
order a 2 fi and the other terms are suppressed by an additional power of a 2 . 
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3.2 Bosonic amplitudes 

In SQED, an electron anti-muon bound state mixes with a selectron anti-smuon bound 
state through photino exchange. To calculate the energy spectrum, there are a number of 
additional scattering diagrams that must be computed. 




= -2ie 2 
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• e±£± e ± /i± 
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(3.14) 
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3.3 Fermionic amplitudes 

In addition to bosonic bound states in SQED, there are fermionic bound states involving an 
electron and anti-smuon or selectron and anti-muon. 
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e p,± — > e-p/x 
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4 Mixing matrices 



In the previous section we have derived the non-relativistic expansions of the potential en- 
ergy V(f,p) between an electron and an anti-muon or their super partners in terms of the 
relative coordinate r*and the relative momentum p. To find the effective quantum mechanical 
description of this system, we also need to expand the kinetic energy of these particles 



£kin = V 7 ™ 2 + P 2 e + V M2 +Pm ~ (™, + M) 



(4.1) 



to the same order, i.e. to fourth order in the momenta. Then the Hamiltonian becomes 



H 



P a 

o + H- mt (r,p) 

2fi r 



(4.2) 



where V and H int are matrices acting on the various "spin" states \s e s m ) where s G {t, 4- 
,+,—}. The components of V are the several potentials given in Sec. [3] In (14. 2 p we have 
singled out the non-relativistic kinetic energy and the Coulomb potential, and denote all 
other terms by H int . 
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We are interested in the bound state spectrum of this system. Without the interactions 
ifint, the solutions would be the familiar hydrogenic bound states \nlmi) with the Bohr 
energies E n = — /ia 2 /2n 2 , see App. [Bj Our task now is to determine the a 4 corrections 
to this spectrum, which have two different sources. The first one is first order degenerate 
perturbation theory. Most of the terms in the scattering amplitudes are of order a 4 and lead 
to mixing between the states in the highly degenerate levels of a given n and I. 

There are a handful of terms in the scattering amplitudes that are of order a 3 , namely 
the first terms in (I3.12p . (I3.13p . (I3.2UI) . and (I3.2ip . Naively, these terms should dominate 
the a 4 contributions. However, it turns out that 



Thus, these terms do not contribute at the level of first order perturbation theory. However, 
as was noted in [5J, they can and do contribute at second order. Recall the formula for the 
energy corrections 



where E\ are the eigen-energies of the bare Hamiltonian. For (i\H- mt \j) of order a 3 and 
E^ of order a 2 the second order corrections will be of order a 4 . The sum in (14.41) should 
be carried over both discrete and continuum states of the hydrogen atom. To carry out the 
sum, we will make use of Schwinger's Coulomb Green's function [9]. 

At first glance, the diagonalization problem of the \nlmi) states seems formidable. For a 
given n, we have n different Z's, for each /, we have 21 + 1 different m^s, and for each mi, we 
have 16 different "spins" |s e s m ) all of the same energy. As it turns out, states of different I 
do not mix. Moreover, the total angular momentum in the z direction is a good quantum 
number. The largest matrix we will need to diagonalize is 6 x 6. 

That the total ^-component of angular momentum is conserved is obvious, but that states 
of different I do not mix is surprising. Both the second order perturbative corrections and 
the hyperfine interaction have the potential to mix an I state with an / + 2 state. For the 
hyperfine interaction, one can check explicitly that the overlap integral 



vanishes. Another integral, which we discuss in Appendix [El guarantees that there is no 
mixing of states with different / in second order perturbation theory. 

In the appendices, we discuss separately the contributions from first order degenerate 
perturbation theory to the / = and I > cases. The reason for the separation is that 
the Dirac delta functions in the scattering amplitudes are only important for I = states 
while the spin-orbit interactions only contribute when / > 0. In the appendices, we also will 
calculate the contribution from second order perturbation theory. Below, we present the 
final result for the mixing matrices. 




(4.3) 




(4.4) 




(4.5) 
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4.1 Overall shift 



The relativistic correction to the kinetic energy as well as a term ~ in the potential V 
do not depend on the spins of the particles. Therefore, these terms lead to an overall shift of 
the levels specified by n and /. We can compute this shift independently from the splitting. 
It is given by the expectation values of the following terms in the \nlmi) basis: 



SE„ 



a 



nl 



Mm \2r 

A 



2 + lf^! + 7r5 (3) (r1 



2r 3 



8\m 3 M 3 ' " 



n 



3 

21 + 1 " 8 + 8Mm 



(4.6) 



4.2 I = 0-states 

There is an eight dimensional space of bosonic bound states with / = 0: 



|in) G {|ft>, lit), III), III), I++), I—), I+-), 1-+)} • 



(4.7) 



(Because the / = sector is already relatively small, we do not take advantage of the fact 
that angular momentum in the z-direction is a good quantum number.) The first entry of 
the state describes the electronic portion of the bound state, whether the electron is spin 
up or down, or whether the selectron comes from the + or — chiral superfield. The second 
entry describes the muonic portion. The Hamiltonian to order a 4 for these states takes the 
form H = E n + 5E n o + Mb- Assembling the contributions from both first and second order 
perturbation theory, the mixing matrix for these states is 



M b = 



2rf 



A B T 



B 




C 




(4.8) 



with 



/ 1 






o 

M 2 +m? 
(M+m) 2 

2Mm 
(M+m) 2 





2Mm 
(M+m) 2 

M 2 +m 2 
(M+m) 2 





\ 





1 



(4.9) 



B 



VMm(M — m) 
(M + m) 2 



-1 
1 



(4.10) 



c 



2Mm 
(M + m) 2 



A curious observation is that M 2 = ^M h . Note that the |H — ) and 
from the other six states; it remains to diagonalize a 6 x 6 matrix. 



(4.11) 



-+) states decouple 
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By supersymmetry, there is also an eight dimensional space of fermionic bound states 
with I = 0: 



in) G |t +), |- t), II +>, h l>, It -), 1+ t), II 1+ I) 



(4.12) 



where the Hamiltonian takes the form H = E n + 5E n Q + Mf + 0(a 5 ). The mixing matrix in 
this case breaks up into a bunch of 2 x 2 blocks: 

M f = ^( D ° 
f 2n 3 V D 

with 



(4.13) 



D 
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M 


y/Mm 








1 






m 








M + m 










M 






\ 








\/Mm 


rn 



(4.14) 



Note, for example, that the |t +) state mixes only with the |— t) state. Again we have 
Mf = ^M f . 



4.3 / > 0-states 

As explained above, there is no mixing between states with different /. Therefore, we fix 
the orbital angular momentum to some / > 0. Furthermore, it is convenient to split this 
space, which contains 8 x {21 + 1) bosonic states and 8 x (2/ + 1) fermionic states, into 
closed subspaces of states with given ^-component, mj, of the total angular momentum. 
The bosonic sector of such a subspace is spanned by the states: 



| in) G {\lmi - l,tt), \lm h tf), |/m ; ,tl), \lm + 1,44), 

\lmi,++), \lmi, — ), \lmt,-\ — ), \lm h — h)| . (4-15) 

There are 21 + 3 such subspaces labeled by rrij — mi — —l — l, —I, ...,/ + 1 where j = I — 1,1, 
or /. States in the set (14. 15j) with magnetic quantum number outside the range —I, / are 
understood to be absent. Thus the dimensions of these subspaces are 1, 7, 8, 8, ... , 8, 7, 1. 

The Hamiltonian acting on these states can be written to order a 4 as H = E n + 5E n i + Mb. 
The mixing matrix takes the form 



21(1 + l)(2l + l)n 3 




(4.16) 



In the leptonic sector it is given by 

/ mi - 1 



M 



A 



■ci 



M 



M+m 



C-'l, 



M+m" l - m ' M+m 
M—m 



M+m 



-mi 



-n i t 



M+m 





mi 







M+m° lm l 



Q,— mi 



M+m 1 
M 



M+m Clm l 



M+m° lm l 
M+m Clm l 

-mi -1 J 



(4.17) 
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and the mixing between leptons and sleptons is given by 



B 




-mi -mi -c imi 



(4.18) 



where Q mi = — mi) (I + mi + 1). There is no interaction among the sleptons, C — 0. 

In the fermionic sector, the states have half-integer total magnetic quantum number 
rrij = mi + | where the range of mi is — I — 1, —I, ... ,1. The corresponding 21 + 2 subspaces 
for j — I ± 1/2 have dimensions 4,8,8,.. .,8, 4 and are spanned by 

| in) e {|Zm { ,t+), \lm h -f), \lm t + 14+), + 1,- I), 
|Zmj,t-), |Zmj, + t), |/m z + l,|-), + 1, + 4,) 



For the mixing matrix in this sector, we find 

4 

iii' 

M f 



21(1 + l)(2l + l)n 3 



D 





D 



with 



/ 



D 



\ 



M+m mi 
sj Mm 



M+m 7711 

M+m ° lm l 

VMm 
M+m 



mi 



C-imi 



VMm 
M+m 

M+m mi 
VMm 



M+m ° lm l 
M+m ° lm l 



M Qj 

M+m Ulm l 

VMm 
M+m Clm l 

M (m, + l) 
(mi + 1) 



VMm 
M+m Clm l 

M+m Clm l 
'Mm 



M+m 

VMm 
M+m 



(mi + 1) 



M+m 



(4.19) 



(4.20) 



(4.21) 



5 Energy splittings and eigenstates 

At order a 2 in the fine structure constant, the energy spectrum is given by the 16n 2 -fold 
degenerate Bohr levels 



-E^nlmiSeSm E'n 



2^2 



(5.1) 



They receive a spin independent shift 5E n i at order a 4 , which we have already computed in 
Sec. 14.11 In this section we calculate the additional splittings of these levels and find the 
energy eigenstates. The splitting energies and the eigenstates are given by the eigenvalues 
and eigenvectors of the mixing matrices Mb and Mf computed above. Because spherically 
symmetric states (I = 0) and asymmetric ones (/ > 0) split up differently into two and three 
levels, respectively (see Fig. [T]on page EJ), we discuss these two cases separately. 

We organize the eigenstates that remain degenerate at order a 4 into multiplets of the 
underlying supersymmetry algebra 



[Jai Jb\ i^abcJc j \JaiQ ! 

[J a ,H] = , [Q a ,H} = 



a, 



a)%Q P 



[Ja,QU = -U<ra)jQL (5-2) 



{Q a ,QU = H5 a , 



(5.3) 
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Figure 2: Action of supercharges. The action of the Q a is indicated by the arrows. Additionally 
there is factor of yfm or \[M when acting on electrons or muons, respectively. The action of Q a is 
the inverse of the action of Q a . 

where J = L + S e + S m is the total angular momentum operator, Q a , a = 1,2, are the 
supercharges, and H is the Hamiltonian. The action of the supercharges on states to zeroth 
order in a is depicted in Fig. |2j To this order they anti-commute to the rest energy m + M. 
We denote super multiplets by IZj where j = 0, ~, 1, . . . refers to the total SU(2) spin of the 
highest submultiplet, i.e. the one whose states are annihilated by the supercharges Q} a . In 
terms of spin-j multiplets Vj of SU(2), the super multiplet IZj is built from Vj_i/2©2Vj©Vj + i/2 
for j > 1/2 and from 2Vq © V\i% for j = 0. The dimension of TZj is (8j + 4). 
The energy eigenstates depend on the mass ratio r = 



5.1 / = 0-states 

The mixing matrices that need to be diagonalized in this case are given in (14. 8[) and ( I4.13P 
for the bosonic and fermionic bound states, respectively. We find that there are 4 bosonic 
and 4 fermionic states with eigenvalue AE = 0, and 4 bosonic and 4 fermionic states with 
eigenvalue AE = It turn out that the former states constitute two super multiplets 7Z , 
while the latter ones fill one TZ\/2- See the / = column of Fig. [1] where these results are 
visualized. 

The energy eigenstates in the first TZq with AE = are given by 

V : |-+) (5.4) 

*V v^ lt+) - \f^Tr\-V , ^r T \l+)- V^M-+> ( 5 - 5 ) 

Vo ■ &(\n) - lit)) + r^(^l++) + I")) (5-6) 
the ones in the second TZq also with AE = are 

Vo : -^(itl) - lit)) + + r\—)) (5.7) 

Vl/2: \f^\t-)-^\+t) , ^ T \i-)-^\+i) (5-8) 

V : |+-) (5.9) 
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and the ones in TZ\/2 with AE = ^3 are 

V* : -&js(\U) ~ lit)) " ?^(l++) " I"")) 



V 1/2 : 



I) 



Vi : |tt) 



1 

V2 



III) + lit) 



III) 



TT7lt-) + A/iT7l+T) , \/tT7II-) + \/tT7I+I) 



(5.10) 
(5.11) 
(5.12) 
(5.13) 



5.2 / > 0-states 



The relevant mixing matrices are (14. 16ft and ( 14.201) . They correspond to the subsector of 
states with fixed principal quantum number n, fixed orbital angular momentum I, and fixed 
z-component of the total angular momentum rrij. As we argued in Sec. HI there is no mixing 
with other subsectors even though sectors with different I are degenerate at order a 2 . 

Super multiplets, however, can only be formed by grouping together states with all possi- 
ble values for rrij. The reason for this is that although a sector with fixed (n, I, rrtj) is closed 
under the action of the Hamiltonian, it is not closed the action of the angular momentum 
operator J nor the supercharges Q a . These latter generators carry spin themselves, and 
therefore can change the m^- value of the state they act on. Thus, we look at all 4 • 4 • (21 + 1) 
states with a given Z-value at once. We find that they form two super multiplets IZi with 



AE = 0, one super multiplet 7^+i/2 with AE 
Hi-x/2 with AE 



2(1+1) (2/+l)n 3 ' 



and one super multiplet 



/10 

~ 2l(2l+l)n 3 ■ 

The states of the first unperturbed IZi are given bj|j 



1 2 



V2I+1 



l + mi + l \mi + 1,1+)- H + 1, - I))] (5-14) 



Vi 
Vi 



\ m h — H) 



j£ (K, tl) - \m h it)) + y+t ( r K> ++) + l m *> — ) 
Jm [v / ^ + ^ + i(v / r^ K t +) - Vr^ K - 1>) 



+ a// - mAJj^ \mi + 1,| +) 



l+r 



\mi + 1,-|) 



(5.15) 
(5.16) 



(5.17) 



3 As discussed in Sec. 14.31 the states are labeled by mi = — I — 1, —I, . . . , I + 1, and kets that end up with 
a magnetic quantum number outside the interval [— I, I] are defined to vanish. 
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The ones in the second are 



1 2 



Vi 
Vi 



V21+1 



y/l-mi (y^ H, t -) - \m h + t>) 

- v^T^TTf,/^ |m, + 1,1 -) - \/X |m, + 1, + (5.1 



|"iJ>H — ) 

(N, tl) - It)) + TT7 (K ++) + r N> )) 

aA+^TT( |m,, t -) - K + t>) 



V2I+1 



+ ^ - m ,^yCZ| mj + 14 -) - \ mi + 1, + |) ) 

4 

The states in TZ[ + i/ 2 which receive a positive energy shift by AE = + 2(i+i){2i+i)n- i are 
V, : 



(5.19) 
(5.20) 

(5.21) 



v /2(i+l)(2Z+l) 



>/(i + mj)(Z-mj + l) H - 1, tt) 



2 M 



1 

v / 2l+T 



- ^Z-mOCZ + mj + l) |mj + 1, U) 

_ (l +mi + l)-{l-m l+ l)r | mjj ^ + (i - mt + i)-(i+ m , + i) T ^ 

-^(K++>-K— >)] 
v / H r ^H r T(v^ K t ±) + K =f t>) 

+ y/l-m^y^ |m, + 1,4. ±) + yX |m, + 1,T|) 



(5.22) 



(5.23) 



- v /2(i+l)(2Z+l) 



^/(Z + m^^ + mj + 1) |m, - 1, tt) 
+ ^/(Z-mOCZ-m^l) |m, + 1, U) 

+ v /(Z + m i + l)(Z-m^ + l)(|m^|t) + H,tl))] (5.24) 
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and the ones in 7^_i/ 2 whose energy is lowered by AE 
Vi-i : 



/ict 



y/2l(2l+l) 



2i(2«+l)n 3 

^ (I - mi )(l - mi + I) \mi - 1, tt) 

+ s/(l + mi)(l + mi + l) \mi + 1, U) 
- y/(l + mi)(l-mi)(\mi,\.t) + K,tl) 



have the form 



(5.25) 



2V i 

' 2 



v / 2l+T 



- ^/Z + m, + lCyiir | m , + 1, 1 ±) + | mj + 1, =f I))] (5.26) 



^21(21+1) 



y/(l+mi)(l-mi + l) \mi - 1, tt) 
- a/(/ - m z )(/ + m ; + 1) |mj + 1, ||) 

+ (i-m,)-g+m,)r ^ _ ( ; +m ,)_(f_ m0T ^ ^ 



+ f^(K++)-K — ) 



(5.27) 



6 Conclusions and outlook 



A comprehensive summary of the results of our computation is given at the end of the 
introduction in Sec. 11.11 Here we discuss some consequences and applications thereof. 

Oscillations Because of the energy splitting, there is an oscillation between different "fla- 
vors". Say we prepare an atom in the flavor state |++) with I = 0, then it can oscillate into 

| ) and :^(|tl) ~~ lit))- The probabilities of finding the atom in one of these states at 

time t after the atom was purely |++) are 



P 



l++> 



P 



I— ) 



\n-w 



1 + 6r 2 + r 4 + 4r(l + r 2 ) cos AE t 

8r 2 (l - cos AEt) 

4t(1 -r) 2 (l - cos AEt) 



(6.1) 
(6.2) 
(6.3) 



where AE = ^p-. If we plug in the actual mass of the electron m ~ 0.5 MeV and the muon 
M « 100 MeV, then AE w 0.7 meV for the lowest state (n = 1). In the limit m < M, the 
oscillations have a frequency of 



u « 10 12 Hz 



m 



500 keV 



137 -«) 4 . 



(6.4) 
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Supersymmetry We can write the mixing matrices in terms of a superpotential W as 
Mf, = W^W and Mf = WW^. The matrix W can easily be constructed from the eigenstates 
given in Sec. as follows. Let A = diag(i?i, . . . , E s ) be the eigenvalues of M b and Mf in some 
fixed order, and let V and U be matrices whose columns are the corresponding eigenvectors. 
Then the superpotential is given by W = U^AV\ If we pair up the eigenvectors in U and 
V appropriately, we can set W = Qi + Q\ or W = Q 2 + Q\- The prerequisite that all E t > 
is satisfied for the I = sector, and in the other sectors we can achieve this requirement by 
adding the identity matrix times the smallest eigenvalue to the mixing matrices. 

Bose condensate Ignoring interactions between these hydrogenic atoms, including any 
instability to form molecules, what happens if we place a large number of these atoms in a 
box? If the atoms are fermionic, then only one fermionic atom can rest in the single particle 
ground state of the box. (More generally, a small but finite number of fermionic atoms can 
rest in the ground state if the ground state has a small but finite degeneracy.) In contrast, 
there is no limit to the number of bosonic atoms that can exist in the single particle ground 
state. Through emission of a photino, a fermionic atom can convert into a bosonic atom. 
Given our assumptions about the absence of interactions, the multiparticle ground state 
will contain at most one fermionic atom. There should be no Fermi sea for these "perfect 
atoms" . It would be interesting to see what changes if any occur to this qualitative picture 
when interactions between the atoms are considered. 

Supersymmetric chemistry For bound states of higher charge nuclei and more than one 
electron, the Pauli exclusion principle will play a much weaker role then it does in traditional 
atomic physics. An electron in an excited orbital can always reduce its interaction energy 
with the nucleus by converting into a selectron and moving into a lower orbital at the possible 
price of increasing its interaction energy with other orbiting selectrons. A supersymmetric 
periodic table should look quite different from the periodic table we are used to. There may 
be additional interesting effects related to these atoms' ability to form molecules. We leave 
a study of such effects for the future. 

Acknowledgments 

We would like to thank Silviu Pufu, Stefan Strieker, Aleksi Vuorinen, and Lian-Tao Wang 
for discussions. We also thank Tomas Rube and Jay Wacker, who independently calculated 
the spectrum of supersymmetric hydrogen |10| . for bringing our attention to a mistake in an 
earlier version of this manuscript. The work of C.H. was supported in part by the US NSF 
under Grant Nos. PHY-0756966 and PHY-0844827. 

A Notation and conventions 

We use the metric rf v = diag(— , +, +, +) and the Levi-Civita symbol e 12 = — ei 2 = 1. 
Implicit contractions of two-dimensional spinor indices are defined as ipx = ^"Xqi $X = 
^aX Q i an d complex conjugation acts as (ip a )^ = ip a , (ip a )^ = ?pa, (ip a XaV = Xa^ a ■ We 
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raise and lower spinor indices from the left: ip a = e^ipp, ip a = e a p^ We employ the Pauli 
matrices 

= , a^ = e^e^a^ = (-l,-a) (A.l) 

to define the Dirac matrices as 

f = ( t < C|' 1 ) • T* - Wf = ( 7 J ) . (A.2) 

The chiral projectors P± = |(1 ± 27 s ) have the matrix representation P + = diag(l,0) and 
P_ = diag(0, 1). The Clifford algebra relations are given by [a li cj v + a u a^) a ^ = — 2r] IMU 5 a 13 , 
(a»a v + o v o»)"p = -2r)» v 5"p, {7^,7^} = -2^1. 

Superfields and components The vector superfield V = —Qa^Q A^(x) + i9 2 9x(x) — 
% 9 2 8x{x) + \9 2 9 2 D(x) contains the photon A^, the gaugino x, an d the auxiliary real scalar 
D. The chiral superfields <P± = <j>±{y) + V2 9ip±(y) + 9 2 F±(y) (where y^ = + i9a^9) 
contain the slepton 0±, the leptons ip±, and the auxiliary complex scalars F±. We introduce 
a Dirac spinor for the leptons 

i> = ( ^ ) , i> = ^7° = ( -r- ) , (A.3) 

and a Majorana spinor for the photino 

A = ( *S ) , A = AV = ( -x° -x Q ) , (A.4) 

subject to the condition A = X c = C*A T with the charge conjugation matrix C = 27°7 2 = 
— C T = — = C* = —C^ 1 = diag(icr 2 , — ia 2 ). The sign in the gauge covariant derivative 
V^X = d^X + iqA^X is determined by the U(l) charge q of X. The field strength is 
F^v = d fl A u — d u A^. 

B Bound states in Coulomb potential 

The free Hamiltonian is given by H = — f. We use the constants \i = for the 

reduced mass, a — ^- for the fine structure constant, and a B = — for the Bohr radius. The 
wave functions ip n i m (r) = (r\nlm) for the bound states are given by 

= ~J^j^ ( —) ' (— ) -p (-— ) *.(«, <p) (b.i) 

a/oI^ 2 y (n + iy. \na B J n L \na B J V na B J 

Their energies are E n = — f^r ■ 111 Mathematica one has to write L b a (x) = LaguerreL[a, b, x] 
and Yi m (9,ip) = SphericalHarmonicY[l, m, 9, ip\. 
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Expectation values 

(nlm\ — \nlm) = , (nlm\ — \nlm) 
r n 2 r 2 



/rcr 



. , ,. , (n00|^|n00) = 
n 3 (Z + \) r 3 



3n 2 



p 2 u 3 o; 3 
(nlm\ — \nlm) = — — 
r rr 



2n 



(nlm\p A \nlm) 



4 4 

/x a 



r?; J 



4n 



/ + 



-3 



(B.2) 
(B.3) 



(nl'm'\ — \nlm) 



/i 3 a 3 
7in 3 



Snn'SloSi^5 m0 S m ' 



ri 6 a 6 



n 3 /(/ + i)(/ + l) 



(nl 



W| (r 1 ^ + 27r(5 (3) (f) | n/m) = 

r 6 rr 



n 



l + \ 



Sll'Sr, 



5 Hi 5 



ll'Omm' 



(B.4) 

for I, V > (B.5) 
(B.6) 



Angular integrals Define components of the unit position vector r as f 
f± = ^ = sinfle^. 



cos 9 and 



, '(I' + m'Ml' -rri) . , l(l + m)(l-m) c 

(Z m \r \lm) = J (V,z+i<-W,m + \\ /r>7/ , ^, nl , ., N <V,i -id m > >m (B.7) 



(2Z' + 1)(2Z + 1) 



(2Z' + 1)(2Z + 1) 



,„ , (Z'±m'-l)(Z'±m') ^ l(frMM A 

(Zm |r±|Zm) = =f^/ — — <ty,/+i<W,m±i ± W , n1 , , ttttt; . -x Oi',i-iO m ',m±i 



(lm\fofo\lm) 



(Z m + l|f + f |Zm) = 



(2Z' + 1)(2Z + 1) 

2Z(Z + 1) - 2m 2 - 1 
(2Z + 3)(2Z- 1) 

(2m + 1)q 

r 



, (lm\r + r-\lm) 



, (Z m + l|f + f + |Z m — 1} = 



(2Z' + 1)(2Z + 1) 

2Z(Z + l) + 2m 2 -2 
(2Z + 3)(2Z- 1) 

2Q m Q 



(2Z + 3)(2Z-1) 
where here and below we use 

ci m = a/(Z - m)(Z + m + 1) . 



(2Z + 3)(2Z- 1) 



(B.8) 
(B.9) 

(0.10) 
(B.H) 



C Feynman Rules 

Fields and particles 



Field 


0+ <t>\ 


<\>- 


0t 








^_ 


X 


X 


creates 


e\ eZ 


eZ 


eZ 
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e + 


el 


el 


7+ 


7- 


annihilates 


eZ e% 


eZ 


eZ 


e + 


e+ 


el 


el 


7- 


7+ 



The U(l) e charges of the particles are encoded in the labels as Q(e q h ) = q, Q(e 9 h ) = q, 
Q{%) = 0, and the U(1) R are given by R(e g h ) = qh, R(e q h ) = 0, R{%) = h. 
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Propagators The momentum p flows from y to x. For photons in Coulomb gauge 
(A^x)A u (y)) A 00 (p) = ± , A 0i (p) = 0, A l3 {p) = --U d - V 4i- 



(ci) 



Coulomb gauge is better suited for taking a non-relativistic limit of the SQED amplitudes. 
For a discussion of the complications that arise when combining Feynman gauge with a 
non-relativistic limit, see for example j8]. For photinos 



(\ a (x)\f)(y)) -> A ap {p) 



+ ie 



a/3 



a/3 



p z — ie 



Electron wave functions 



d 3 p 1 



(2tt)3 ^2Ep 

Contractions with external particles are given by 

(e-(p,s)\i; = u s (p) e~ ipx 
(e + {p, s)\i> = v s (p) e~ ipx 

It is useful to know the identities 

cr • p + m 



btu s {p)e ipx + d s Jv s (p)e 



j—ipx 



V2W+ 



y/o^p 



1>\e-(jp,s)) =u s {p) e ipx 
$\e + (p, s)) =v s {p) e ipx 

a ■ p + m 



m) 



v/2(j9° + m) 



(C.2) 



(C.3) 

(C.4) 
(C.5) 

(C.6) 



In the quantum mechanical setting we convert the spinor 77 for an anti-particle into a spinor 
for a particle using the relation 77 = C£* where the charge conjugation matrix is C = io 2 . 



D Degenerate Perturbation Theory 
D.l I = 

The only non-zero matrix elements of the interaction Hamiltonian 

f4 / j 1 



V 



Hdeg — H int + 

for s-wave states |n00 s e s m ) are 

(n00±±|# deg |n00 £ ei £ mi ) = =F 



m 3 M 3 J 47r Mmr 



3 >mo>eo 



2(m + m) + 3Mm at ' a '" 



1 /i 3 a 4 M — m 

2 (.U/,/)^ v s " ' 



1 LL Oi M — 171 
(n00UUo\H dcg \n00±±) = ±- ^ 3 (Mm)3/2 4^C 



(r700±±|i/ deg |n00±±) 



/i 3 a 4 2 
77 3 Mm 
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(D.l) 

(D.3) 
(D.4) 
(D.5) 



for bosonic atoms, and 



(n00£ eo ±|#dcg|™00£ e ;±> 

mi) 



(n00+£ TOO |F deg |n00£ ei ±) 



mi) 



1 



Mm 2m 2 



Mm 2M 2 



l/Za 4 M + m 
2 n 3 (Mm) 3 / 2 

1 /i 3 a 4 M + m , 

2 n 3 (Mm) 3 / 2 ^™ 



for fermionic ones. 



(D.6) 
(D.7) 
(D.8) 
(D.9) 



D.2 / > 

States of different n and / do not mix, and we fix n and /. We begin with the bosonic states 
(I4.15p . The only non-zero matrix elements of the interaction Hamiltonian are 



(nlm'i ieoimo\H Aeg \nlmi £ e iimi) 



— ( In/, 



2M 2 Mm 



1 



{nlm'i ±±\H deg \nlmi £ ei £ m ;) 
(nlm'i ieoimo\H deg \nlmi ++) 



2Mm 
C M + m 



(3f • a e r ■ a m - a e ■ a 
L 



lm, 



T 4 (Mm) 3 / 2 {lm ' 1 



where 



C 



C M + m i 
: T(Mm) 3 / 2 \ 



/i 3 a 4 
n 3 Z(Z + ±)(Z + 1) 



lm. 



Imi) ■ (Cniv 2 °iei) , (D.ll) 
Im^-^lai^Cmo) , (D.12) 

(D.13) 



The nonzero matrix elements of the interaction Hamiltonian for the fermionic states ( 14.1 9p 
are 



(nlm'i i eo +\H dcg \nlmi £ ei ±) 
(nlm'i ±Cmo I H deg | nlmi ±£ mi ) 
(nlm^ T€mo\Hde g \nlmi £ ei ±) 
(nlm[ ieo+\H Aeg \nlmi ±£ mi } 



C 



1 1 

+ 



2 VMm 2m 2 

C / I 1 

+ 



2 \Mm 2M 2 
C M + m 
T (Mm) 3 / 2 
C M + m 



(Zm;£ eo ±|£-c?|Zm z £ei±) , (D.14) 
(/m^±£ mo |L-a|/m;±£ ?ni ) , (D.15) 
(Jm{|£|Zm,).£L*&< , (D.16) 



4 (Mm) 3 / 2 
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(Zm{|L|Zmz) ■ £[ a£„ 



(D.17) 



Thus to determine the energy corrections from first order degenerate perturbation theory, 



we need the matrix elements of L ■ a, a e ■ a m , and f • a e r ■ a m : 







(L-a) 


|m/ t) |mj + 1 |) 








(mi t 1 










(mi + 1 || 


Cim, -(mj + 1) 






(B e • cr m ) 


\mi - 1 tt) \mi It) |m z tl) 


i i 1 i i \ 
\mi + 1 II) 




1 1 A A I 

(mi - 1 tt I 


1 










/ i^l 

(m It 1 





-1 2 


o 




\' lb l 1 4- 1 





2 -1 


n 




(m, + 1 || | 








1 


(r ■ a e f ■ a m ) 


|™J - 1 tt) 


\mi |t) \mi tl) 


\mi + 1 ||) 


{mi - 1 tt 1 


cos 2 9 


sin ^ cos 9 e~ l< ^ sin ^ cos 9 e~' 1 ^ 


sin 2 9 e-' Mct> 


(m It 1 


sin 9 cos # e l 


^ — cos 2 


9 sin 2 9 


— sin 9 cos # e~ 1 ^ 


(m tl 1 


sin 9 cos # e l 


* sin 2 6» - cos 2 9 


— sin 9 cos 5 e~ 1 ^ 


(m t + 1 || | 


sin 2 9 e 2i(t> 


— sin # cos 9 e 1 ^ — sin 9 cos # e w 


b cos 2 


evaluating this last matrix, the integrals (IB.9HB.10 ) are useful. We also find 




±{ia 2 L ■ a) 


l^i - 1 tt) 


\m t |t) |mj tl) 


771; + 1 ||) 




(m + +| 




— mi —mi 


— Qm; 




(mi 1 


Q,— m; 


mi mi 





(D.18) 



(D.19) 



(D.20) 



(D.21) 



E Second Order Perturbation Theory 

We are interested in computing the correction to the energy of a state \nlm, s e s m ) at second 
order in perturbation theory. (For ease of notation, we remove the subscript I from mi in 
this subsection and replace m with m e .) These second order corrections will not mix states 
of different n and I, n because the energies are different and I because of the vanishing of an 
integral we discuss below. We need to compute the matrix 

(nlm, s e s m \H int \i)(i\H int \nlm', s' e s' m ) 



AE(n, I, m, mf, s e , s' e , s m , s' m ) 



E„ — Ei 



(E.l) 



where the ' on the sum means we should omit the states with Ei = E n . (We can omit these 
states in the sum because of eq. (14. 3p .) This sum involves both the discrete and continuum 
hydrogenic states. 

The Hermitian matrix AE has a block diagonal form. The four bosonic states 



\m - 1 tt), \m It), \m tl), \m + 1 ||) 
mix among themselves according to 

/ AE_+(m-l) -AE_ (m) 

— AF,„,(rr,.— -\\ 



(E.2) 



AE 0+ {m-i) 
AE 0+ {m-l) 
\ -AE ++ (m-l] 



AE 00 (m) 
AE 00 {m) 
AE +0 {m) 



AE 00 {m) 
AE 00 {m) 
AE +0 {m) 



-AE_ 
AE . 
AE . 
AE, 



.(m + 1) 
m + 1) 
m + 1) 
(m+1) J 



(E.3) 
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The factor of two comes from the sum over scalar intermediate states ++ and . We will 

define AEij(m) presently. The coefficients of AE vanish for the other four bosonic states 

\m, ++), | m, — ), |m, H — ), | m, — h) , (E.4) 

provided I > 0. In the special case I = 0, the states \m, ++), \m, ) mix as 

AE r AE r 



AE r AE r 



(E.5) 



where AE r = AE 00 + (AE + _ + AE_ + )/2. For the eight fermionic states, AE reduces to 
four 2x2 blocks. The states \m, f ±), \m + 1, \. ±) and the states |m, ± f), (m + l, ± \) each 
have the same second order mixing matrix 



AE 00 (m) + AE_ + (m) AE_ (m + 1) - AE _(m + 1) 
AE 0+ (m) - AE +0 (m) AE 00 (m + 1) + AE+_(m + 1) 



(E.6) 



The orbital part of the interaction that contributes at second order in perturbation theory 
has the schematic form (p — q)i/(p — q) 2 . We will work in a basis where po = p z and 
p± = (p x ±ip y ). 

To compute the sum (lE.lj) . we make use of a beautiful result of Schwinger [9] for the 
Coulomb Green's function in momentum space: 

(e - Q G(p,p>; E) + ^J J 0^G(P",P'; E) = f) . (E.7) 

Introducing a p such that E = —pl/2fi, the construction makes use of a mapping of the 
four momentum (po,p) to an S 3 : 

£=-^Sv, &> = 4^> r 2 + e 2 = i - (e.8) 

P +Po Po+P 

Let f2 = (-0, 6 1 , 0) be a set of angular coordinates on the S 3 . The matrix elements in the 
Hamiltonian at second order in perturbation theory can be built from the expression: 

at? f \ 47r2ft2 f d 3 pd 3 p'd 3 qd 3 q' (g - p'), , , (q - q') j 

(E.9) 

where 

n(v J J, t?\ 16 ^PO \^ Y n'l'm'{^p')yn'Vm' < S l q') (T?in\ 

(j{P A,E) = — > — . (h.IUJ 

[Pt+P 2 r(Po + (l ) ^ l-afi/n'po 

The conserved energy is E — —fia 2 /2n 2 . The <p n i m are the hydrogen wave functions in 
momentum space: 

4(>o) 5/2 

4>nlm{p) = 7^5— 2T?^n/m(^) ! Y n i m (Q) = Z ni (lp)Yi m (9, if) . (E.H) 

[P 2 + P0) 2 
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The Y n i m (Q) are spherical harmonics on the momentum S 3 . The Y n im and Y\ m satisfy similar 
orthonormality and completeness relations. More explicitly, 



ZM) = N nl sin^C^.^cos^) 



where the C l n (x) are the Gegenbauer polynomials and 



N, 



nl 



n(n-l-l)\2 2l+1 {l\f 



{n + l)\ 



71 



1/2 



(E.12) 



(E.13) 



Making use of the relation 

1 1 



4tt 2 (£ - e) 2 



(E.14) 



nlm 



the expression (1E.9I) simplifies to 



/i 3 a 4 1 

Mm„ n 2 ^— ' n' — n~'' 



" Qn'l'rn 1 ,m+i+j Qti'l'm' ,m* 



n'l'm' 



where 



5, 



n'l'm' m " " 



>' — n f £■ 



To perform this integral, we use the definition of the F n ; m (f2): 



5 



n'l'm' m n n 



2nn' 



sin t/> 



r (Zm|ri|Z'm') / -Z n , v (i})) dip . 

' J 1 + cos 



(E.15) 



(E.16) 



(E.17) 



The integral (/mlfjl/'m') vanishes unless / = /' ± 1 by the electric dipole selection rules. The 
energy shift simplifies to 



,3^,4 



/j roc 
Mm e n 4 



n'V 



n — n 
An' 



1 + cos 1p 



where c± depends on i, j, I, and m but not on n', I', or m! . 
Let's define 

I{n,n',l,l') = I" ^^T c'n' +V-i(cos^)C^_i(cos^) # • 



/ l + cos^ 
For positive integers n and n', we find that 



(E.18) 



(E.19) 



\I(n, n', I, I + 



(f n(n+l)\ \ 2 i 

1)| 2 = { ^Wi!(i + i)!(„-!-i)! j ' n < n 



n > n' . 



(E.20) 
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Given these results, we can evaluate the sum for / > 0: 



AE i3 (m) = 



/j, 3 a 4 



Mm e n A 



Mm P n 3 



n(n — I — 1)! 
(n + l)\ 

n(n + l)\ 



2>~ w > (n'-l)l 



n'=l 

oo 



n-l-l)\ ■ 

' n'=n+l 



— c_- 



^ [n n) (n' + l + l)\ 



(2Z + l)(2Z + 2) "2/(2/ + l) 



(E.21) 



(If / = 0, the coefficient of c_ will vanish because the state V = 1 — 1 does not exist.) 

From the integral /, we deduce that there is no mixing between / and / + 2 states at 
second order in perturbation theory. Note that I(n,n' ,l,V) = I(n',n,l',l). To get mixing 
between these states we need some amplitude to scatter from an nl state to an n', I + 1 state 
and back to a n, I + 2 state. In other words, the product I(n, n' , I, I + l)I(n', n, I + 1, 1 + 2) 
should not vanish. However, I(n', n, 1 + 1,1 + 2) will vanish unless n! < n while I(n, n', 1,1 + 1) 
will vanish unless n < n' . 

To evaluate the AEij completely, we need 



(Imlfill'm^il'm'lfillm) 

i'=-V 
V 

(lm\z\l'm')(l'm'\z\lm) ■ 



2/ + 1 
1 



m'=-V 
V 



+ (i + 1) , 

{l + lf-m 2 



l 2 -m 2 
21 - 1 



H'l-l 



+ 



21 + 3 



i+i 



(E.22) 
,(E.23) 



(l,m±l\r±\l'm'}(l'm'\r T \l,m±l} = 



m'=-V 



(l + m)(l + m + l) (Z + m)(l + m + 1) 

-<V,J-1 H /ot ; lWo7 ; oN — 



(2Z + 1)(2Z-1) ' (2/ + l)(2/ + 3) 

')(ZW|f T |Z,m± 1) = ^ (l,m±l\r±\l'm')(l' m'\z\lm) 

I ± m 



(E.24) 



m'=-l 



Q,±m 
2Z + 1 



=F 



2/ - 1 



m'=-V 

I + m + 1 
21 + 3 



(Z,m ± l|f±|ZW)(ZW|f±|Z,m =F 1} = 



QtnQ,-m 

' 2Z + 1 



2Z - 1 2Z + 3 



(E.25) 
(E.26) 

(E.27) 

AE T ±(m) = (1(1 + 1) - 3m 2 ± (21 - 1)(2Z + 3)m)N , (E.28) 



m'=-V 

The coefficients of Syj-i and Syj+i are c_ and c + respectively. Finally, we find 

AE 00 (m) = (-1(1 + 1) + 3m 2 )N , 



AE 0± (m + 1) = AE T0 (m) = (+21(1 + 1) + 3m)c ltTm N , 
AE ++ (m - 1) = AE__(m + 1) = 3c, m c, _ m JV , 



(E.29) 
(E.30) 
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where 



/x 3 a 4 1 
= ~2Mm e n 3 1(1 + 1)(2Z - 1)(2Z + 1)(2Z + 3) ' ^ 



In the special case I = 0, we find that 



^ ±T (0) = 24^(0) = -J™\ (E.32) 



while Z\£ ±0 = AE ± = 0. 
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